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Abstract 

In this work the semiclassical model of pure Coulomb excitation was applied to the breakup of 
15,17,19q_ 'pj^g ground state wave functions were calculated in the particle-rotor model including 
core excitation. The importance of interference terms in the dipole strength arising after including 
core degrees of freedom is analyzed for each isotope. It is shown that Coulomb interference effects 
are important for the case of ^''C. 
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I. INTRODUCTION 



The study of neutron rich hght exotic nuclei has been the subject of extensive experi- 
mental and theoretical work for more than two decades already. One of the tools to probe 
the structure of these nuclei are the experiments on the Coulomb excitation and breakup 
whose reaction mechanism is very well known. 

Several experiments have been recently performed to study carbon isotopes i5,i7,i9q^ 
was found m that is a candidate to have one-neutron halo. Ground state with spin-parity 
r fiaal „.e-neut.on separation e„e.,v less ... 1 MeV QQS favo. tKe fonnat.n 
of the halo in this nucleus, which is also confirmed by narrow momentum distributions of 
i^C following ^^C breakup [ij. 

There were several experiments to measure separation energy of ^^C. The experiments 
using time-of-flight techniques suggest small separation energy, that is, weighted average 
yields 242±95 keV iJji]. The Coulomb dissociation of ^^C was studied by Nakamura in 
The analysis of angular distributions of breakup products suggests the value 0.53±0.13 MeV. 
Using this value in the simple cluster model calculation of the dipole strength gives good 
agreement with the data. But the analysis of recent experiment of Maddalena et al. 0] 
on nuclear breakup of ^^C yields 0.65±0.15 MeV and 0.8±0.3 MeV. The adopted value of 
one-neutron separation energy of ^^C is given in Ref. yiand Sn = 0.58 MeV. 

Unlike '»C the =pin-parity ot "C was fonnd to be f + fl, with binding energy also smaller 
then 1 MeV p]. However the spin | is an indication of a (i-wave single particle structure 
and accordingly the existence of centrifugal barrier in this case does not favor the halo. 

^^C again has a ground state 0| with a large amount of s-wave in the wave function 
with the separation energy of around 1.2 MeV j^. Although this isotope is not generally 
recognized as a halo-nucleus, there is an evidence that this nucleus has a halo-like structure 
(see 13] and 0). 

Coulomb breakup of ^^C and ^^C at relativistic energies has recently been studied in Ref. 
0]. The main ground state configuration of ^^C is found to be ^■*C(0^)(S)z/s with spectroscopic 
factor consistent with earlier studies. The predominant ground state configuration is of ^^C 
is found to be ^^C(2+) ® Us^. 

A number of theoretical studies of Coulomb and nuclear breakup of heavy carbon isoto pes 



were performed in the framework of post-form distorted wave Born approximation llO, lU 



1 isoto pes 
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12 1 . The study of Ref. la confirmed tlie existence of one-neutron halo in ^^Be/^C and ^^C 
but not in ^^C. In Ref. Coulomb breakup of ^^Be and ^^C on ^osp]^ ^^^g analyzed, where 
the transitions to excited states of the projectiles core were calculated. It was found that the 
contributions of the core excited states are very small. In Ref. the respective roles of the 
first order and higher order effects for different beam energies were investigated. It was found 
that higher order effects are small in case of higher beam energies and forward scattering, but 
important at incident energies < 30 MeV/nucleon. In Ref. a semiclassical dynamical 
model of projectile excitation was applied to nuclear and Coulomb breakup of ^^Be and ^^C 
on ^^'^Pb. It was shown that Coulomb breakup dominates relative energy spectra around the 
peak region while the nuclear breakup is important at higher relative energy. The higher 
order effects are found to be generally small and dependent on the theoretical model. 

In this paper we use the semiclassical model of pure Coulomb breakup. In this model the 
Coulomb breakup cross-section is calculated using the virtual photon numbers formalism and 
the dipole strength function dB{El)/dE which requires the knowledge of the ground state 
and excited states wave functions. For the ground state we use wave functions obtained 
in the particle-rotor model with core excitation. For the continuum states we use plane 
waves. The inclusion of core degrees of freedom gives rise to the interference terms in the 
dipole strength function. The model employed here was used before in Ref. Q], where 
the relative energy spectra were calculated for the Coulomb breakup of ^^C and integrated 
cross-sections for the Coulomb breakup of i5,i7,i9q -^yg^-g presented for different ground state 
scenarios of these carbon isotopes. The lack of experimental data made it difficult to draw 
final conclusions about structure of these heavy carbon isotopes. 

The purpose of this paper is to give a detailed analysis of the particle-core result. In 
particular we assess the importance of the interference terms in the dipole strength function, 
linear in the core deformation, on the shape and value of the cross section. In our approach we 
shall analyze the nuclear-corrected data, namely the ones where the nuclear contribution has 
been removed through the commonly used extrapolation procedure. For a recent overview 
of this scaling method see Ref. 
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II. CORE-PARTICLE MODEL WITH CORE EXCITATION 



The inclusion of core degrees of freedom as it is done in the core-particle model with core 
excitation is simple and physically transparent. This method was described in Q]. It was 
18, [3] to study the properties of such weakly-bound and unbound two- 



employed in 



17 



body systems as ^^Be, ^^C, ^'^Li and three-body ^^Be. The inclusion of core excitation made 
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it possible to describe Be as an s-intruder nucleus. Esbensen et al. [20[ also applied this 



model to positive parity states in ^^Be and ^^C. In the work by Ridikas et al. 1J| different 
scenarios were proposed for the g.s. structure of i5,i7,i9q Recently this model was applied 
io the simultaneous description of borromean nucleus ^^Be and its binary constituent ^^Be 
2ll |. ^^Be was also studied in the particle- vibrator model by Vinh Mau in [2^ . 

The basic idea of this method is that deformation of the core can lead to couplings with 
excited states of the core. The total Hamiltonian H of the two-body system(core -|- n) can 
be written in the following way: 



H = f + Hrot + V, (1) 

where T is kinetic energy of the relative motion of core and valence neutron, Hrot is the 
Hamiltonian of a deformed axially symmetric rigid rotor, V is the interaction between the 
core and the neutron. The wave function of the system with total spin J has the following 
form: 

Ijl mimgmjmj 

{jm,Imj\JM)Yim,if)X,mS^)<l)i^AO, (2) 

where Xijiij) ^i-re the radial wave functions, Yimiif ) are the spherical harmonics, Xsms(^) are 
the spin functions. The core states 0^/0 (0 a-^e eigenvalues of the Hrot and are proportional 
to the rotational matrices: 




where ^ are the Euler angles, characterizing the orientation of the core in the laboratory 
system and I = a/ (2J -|- 1). 

The interaction term in the Hamiltonian consists of the deformed Woods-Saxon potential 
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and standard undeformed spin-orbit part: 



T AUIS 

V-^ir,e')= (3) 
1 + exp — 



V--(r) = 2(— )^-^— a+exp . (4) 

m,rC r ar I a ) 

We include only quadrupole term in the expansion of the radius: 

R{9') = Ro^l + PY2o{cos9')'^, (5) 

where /3 is the quadrupole deformation parameter and 6' is the polar angle of the particle 
in the body-fixed coordinate system. 

To solve the problem we substitute the expansion (2) into the total Schrodinger equation 
to obtain the set of coupled equations: 

(T + \A;(r) -E + ej) x^(r) = - ^ V;^^, (r)x^, (r), (6) 

where 7 = {l,j,I}, kinetic energy T = —h^/2fi [d? /dr'^ — l{l + 1)/t^] and matrix elements 
V^y (r) = (/, j, J; JM\V{r, 9')\l',j', /'; JM). ej is the excitation energy of the core state with 
spin /. 

The deformed part of the interaction is then expanded in terms of the Legendre poly- 
nomials in order to separate the angular and radial part in the calculation of the matrix 
elements: 

V(r,ff) = J2VQ(r)PQ{coBff) (7) 

n 

The corresponding angular part of the matrix elements could be found in |16l |. 

The alternative and further simplified way to calculate matrix elements is to expand the 
deformed Woods-Saxon potential in Taylor series over deformation parameter j3 keeping only 
the leading term: 

V{r,e') = Vof{r) + (3v{r)Y2o{cosd') (8) 

/(O = , ^ ^ r-B, (9) 

1 + exp 



where v{r) = Vo/(r)^exp {^—^) ^ is the radial part of the coupling matrix element. This 
method allows the analytical calculation of the radial part of the coupling matrix elements. 



TABLE I: Angular matrix elements cai for the ground state ^ (^^C/^C) and | (^''C) 



1 + 

2 


3 + 
2 


i 


1 


2 


3 


i 


1 


2 


3 


1 





0.218 


0.178 


1 





-0.126 


-0.126 


2 


0.218 


0.14 


0.044 


2 


-0.126 





0.126 


3 


0.178 


0.044 


0.12 


3 


0.126 


0.126 






In our analysis we will employ particle-core method with core excitation to describe ^^C 
and ^^C. We also use this model to describe ^^C, although the core ^^C is not a rotor. It 
is already established that the ground state of the ^^C is 1/2+ state. Thus, for the ground 
state of ^^C and ^^C we will have to couple following three channels: 
1 + 

I- ) = ai|2si/2 ® 0+) + a2|14/2 ® 2+) + Ogll^/a ® 2+), 

For 1^ ground state in ^''C we also couple three channels for simplicity: 

3 + 

I- ) = fei|14/2 ® 0+) + 62I14/2 ® 2+) + 63I2S1/2 ® 2+). 

The system of coupled channel equation using potential expansion of Eq.® after calcu- 
lation of angular matrix elements (/jJ; JM\Y2o{cos9')\l' jT] JM) has the following form: 



(T, + Vr{r) + V,f{r) + Cu(3v{r) + - E) x.(r) = 

-(3vir)J2cu'xAr), (10) 

where index i numbers the definite channel \ljl), Vf°{r) is a spin orbit potential in channel 
i and Cu' is array of angular matrix elements (see Table I!]). 

The solution of the system of coupled channel equations Eq.(jni) or Eq. (|Tnjl gives the bound 
state energies and the corresponding radial functions of the components with their relative 
stren gth s. In the present calculation we used the R-matrix method on Lagrange mesh (see 



Ref. 



23I ] and references therein) to solve the system of coupled channel equations 
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III. COULOMB BREAKUP CROSS SECTIONS AND DIPOLE STRENGTH 
FUNCTION 



In the semiclassical model of Coulomb excitation and breakup the process of absorption of 
radiation is treated quan tum mechanically but the projectile is assumed to move in a straight 
line (see Refs. jsJj ^^^^ order perturbation theory the process of Coulomb 

excitation can be described as emission and absorption of virtual photons. Thus, Coulomb 
dissociation spectrum is related to the dipole strength function dB{El)/dE through 

ddc 16n^dB{El) 1 



POO 

/ 2nbdbNEi{E*,b), (11) 

Jbo 



dE mc dE E 

where E is the relative energy between the core and the neutron, E* is an excitation energy, 
NEi{oj,h) is the number of virtual photons, u = E*/h. bo is the cutoff parameter, which 
is approximated by the sum of projectile and target radii. Here we assumed sharp cutoff 
theory for Coulomb breakup. Within the framework of a direct breakup mechanism, dipole 
strength function is given by the transition matrix element 



dB{El) 



dE ^ 

M 



I - " Ze 

|<qM|<l>(r)> , rf=— rFiM, (12) 



where < q| represents the scattering state, ^'(r) is the ground state wave function of the 
nucleus and d is the electric dipole operator. 

The cluster model or Yukawa+plane wave approximation j3] is usually used to calculate 
dB{El)/dE distributions. In this approximation one assumes Yukawa wave function for the 
ground state and a plane wave for a continuum state and dipole strength function takes the 
form 

"Bm-SN.^i^Xf-f'l (13) 



dE "7rVV^c + l/ (E + Eb)^ 

where 5* is the spectroscopic factor for the 2si/2 state and A^o is the normalization factor 

see Refs. 0,HQ)- 

Eb is the binding energy of the neutron in the ground state, is the 
reduced mass of core+neutron and Ac is the mass of the core. It is seen that the the shape 
of the distribution depends only on the binding energy Eb and has a peak at E = ^Eb- 
The integral, B{E1) is given by 

The validity of the Yukawa+plane waves approximation is based on the fact that for low 
excitation energy dB[El)/dE is determined by the outer part of the wave function. 



When allowing the core to be excited (deformed or vibrational), then the cross section can 
be expressed as the incoherent sum of components da{I)/dE corresponding to different core 
states with spin I populated after one neutron removal. Furthermore for each core state the 
cross section is further decomposed into an incoherent sum over contributions from different 
angular momenta j of the valence neutron in it's initial state. Accordingly one has the 
general expression for a final plane wave (n+core), (q|, 

da{I) 



dE 

E 



^NEi(E)Y,C^S(I,nl]) 

j 

) I\^rYim\^pnij{r) ® /) 



X 



(15) 



where the value of q depends on the energy of the core state. 

The above form for the exclusive Coulomb dissociation has an analogous, one-neutron 



3G 



removal cross section (nuclear) obtained within the eikonal approximation by Tostevin 
3l| and used extensively for spectroscopic study of radioactive nuclei. 

When the inclusive Coulomb dissociation is required one, namely the final state of the 
core is not observed, then the cross-section becomes 



da 
dE 



Ed(y{I) 
dE 



(16) 



The question we ask here is how large are the interference terms present in the individual 
da{I)/dE and in the sum, Eq (fT^ ? These genuine quantal interference terms arise natu- 
rally when calculating the matrix elements {q\^rYim\ipnij) using the well known Bauer's 
expansion of the plain wave into partial waves 

4-7r 



(q|r) 



kr 



(17) 



Lm 



For each core state / the dipole strength wave function is calculated using the wave 

functions of Eq.Q for the ground states (with spin J) and plane waves for continuum states 

I I 

(with spin J') and has the following form jl7i] : 



dB{El,I,J ^ J') _ P_ 
dE ~ 



5Z ■^ki'f'^ki'j' 



(18) 
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where angular part of the matrix element is given by: 

ZcC I 3 




J I J- » 

J 1 1 / l' 1 1 \ 

(19) 



and the radial part of the matrix element: 



Ki'f = y^l xUr)T3v{kT)dT, (20) 

where Xijiij) ^"^^ g-S- radial wave functions of Eq.Q and jvikr) are the spherical Bessel 
functions. 

The expression for the dipole strength function of Eq. ()18|) contains the coherent sum of 
contributions from different channels (but with the same core state), which appear because 
of the expansion of the total wave function of the ground state in terms of the core states. 
This gives rise to the interference terms in the dipole strength and consequently in the cross 
section. 

To compare with experimental data, one needs to sum the expression p8|) over allowed 
final angular momentum J' . The spin of the ground state of ^^C and ^^C is J = |^ and 
for the final state we have J' = {| , | }. For the g.s. of ^''C with J = \^ one also has 

~ {2 '2 J"' 

The contributions of different terms entering Eq. pHj) will be analyzed in the next Section. 

IV. APPLICATION 

In this section we present our calculations for the Coulomb dissociation of the carbon 
isotopes ^^C, ^^C and ^^C in the field of ^°^Pb. The parameters of particle-core model used 
here for each of the carbon isotopes are given in Table ITTl 
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The role of the transitions to the core excited states in the Coulomb breakup of ^^C (and 

, where the integrated partial cross sections 
and momentum distributions for the ground state as well as excited bound states of core 



TABLE II: Parameters of the model used to describe C, C and C and the corresponding 
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ground state wave functions. 
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depth Vq° and deformation parameter /? are 


the same as m Ket. 14J. 
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V n \ 1V±C V J 


49 


59.33 


71.12 


R 
P 
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1.2 
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7.1 


Sn (MeV) 


0.65 
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1.22 


results 








g.s. structure 


1/2+ 


3/2+ 


1/2+ 


(%) 


71.5 Si/2 0+ 


14.5 d-s/2 <S) 0+ 


87.2 si/2 0+ 




25.3 d5/2 ® 2+ 


76.5 Si/2 2+ 


10.9 (i5/2 (E> 2+ 




3.2 ^3/2 2+ 


9.0 d3/2 ® 2+ 


1.9 4/2 2+ 



nuclei were calculated within the finite-range DWBA as well as within the adiabatic model 
of the Coulomb breakup. It was found that the transitions to excited states of the core are 
quite weak and the interference effects are suppressed. 

In our calculation, there are 7 coherent contributions to the — > | transition and 
4 coherent contributions to the — ^ | transition. In fact, in the |^ — > | case, we 
found out that the dominant contribution is the direct Si/2 — > P3/2 one (see Fig^). All 
the other 6 contributions add up to a 6% of the total, confirming the result of Ref. jlol |. 
The contributions of core excited state 2+ are very small and in order to make them visible 
they are multiplied by 100 in FigH^ and FigUlD and by 10 in FigHt. In the |^ — > | 
transition, besides the dominant si/2 — > pi/2 one, there is only one interference term which 
add up destructively, namely 2(1^5/2, 2+|(i|p3/2, 2+; fc) (1^3/2, 2+|(i|p3/2, 2+; fc). The other 2 
contributions, add up to about few percents, as Fig|TjD shows. Total contributions of all 
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FIG. 1: Contributions to the dipole strength function for ^^C. (a)Dominant contributions to the 
dipole strength for — > | transition; (b) All contributions to the dipole strength for |^ — > ^ ; 
(c) Comparison of overall contribution of squared terms(contribution of 0^ is not shown here) and 
interference terms; (d)overall comparison of contribution of 



I and and | 



^ and 



their sum for 0^ contributions (thin lines)and 2^ contributions (thick lines); 



squared terms (core excited components only) and interference terms are compared in FigO^. 

The total — > | and — > | transitions are shown in Fig^ together with their 
sum, to be compared with the data. 

When compared to the data, the major low-energy peak is quite nicely reproduced, with 
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FIG. 2: ^^C+^'^^Pb at 67A MeV Coulomb breakup cross sections calculated in different models. 
Solid line stands for the cluster model calculation with 5=0.67; dashed line - present calculation 
with £2+ =1.62 MeV; dotted line - present calculation with £2+ =1.2 MeV; dot-dashed line - toy 
model calculation (see text). Data from j^. 

70% s-contribution and 30% d-contribution. The data also exhibit some structure at higher 
excitation energies, which can not be accounted for by the particle-rotor model. In fact, 
we have already taken the liberty of reducing the energy of the 2+ state of the core (since 
the rotor is found inside the halo nucleus and changes may ensue) to see whether this, with 
accompanying change in the relative energy of the halo neutron, can reproduce better the 
data. We found out that the combination €2+ = 1.2 MeV and Sn = 0.65 MeV gives the best 
account, but still misses the strength at higher energies. Since the total contribution to the 
cross section involves the incoherent sum of the |^ — > | and — > | , and the former 
exhibits a significant structure besides the main peak, we were tempted to consider a two- 
cluster model for the cross- sect ion, one peaks at E = ^Eb, while the other at E = ^{X—Eb), 
with X related to the excitation energy of the 2+ state in the deformed core. The result of 
our calculation is shown in FigEl The "spectroscopic factors" attached to the two cluster 
response are, respectively, 0.7 and 0.3. The above toy model is instructive, as it points to 
ways of improving the particle-excited core model. Of course, the Coulomb dissociation, 
being, what is known as, elastic breakup (no target excitation) may have to be corrected 



owing to possible contributions of nuclear- Coulomb interference 13| and inelastic breakup 
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B. 



In the dissociation of ^^C, whose ground state spin is we again sum incoherently two 
transitions |^ — > | and |^ — > | . Unhke the case of breakup of ^^C, here we have ^'^C 
mainly with the core excited contribution in the ground state wave function (see Table |n)) . 
Therefore terms c/3/2 — > P3/2 and ^3/2 — > pi/2, which corresponds to the contribution of 
the 0"*" in the core, are now small. 

For 1^ — * I there are all together 10 contributions of which 6 are interference 
terms (see FigOK)- For I"*" — > | transition, the dominant terms I"*" — > \ and 
2(2si/2, 2~''|(i|p3/2, 2+; /i;)(l(i3/2, 2+|(i|pi/2, 2+; fc) cancel each other. As FigOlshows, the over- 
all strength is predominantly due to the 2"*" state in the core and is rather small because 
of the cancellation of the two dominant terms above. The contribution of the O"*" in the 
core, namely | (1(^3/2, 0"'"|(i|p3/2, 0+; A;) p is but a few percent. This shows clearly that in the 
I"*" — > I transition quantum interference and core excitation are very important. 

The transition — > ^ contains altogether three squared and one interference terms, 
the dominant squared one is |(2si/2,2+|c?|p3/2,2+; and the only one interference term 
2(2si/2, 2+|(i|p3/2, 2+; /c)(l(i3/2, 2+|(i|p3/2, 2+; /c) positive and almost equal the direct term 

•^3/2 — ^ Pl/2- 

Transitions |^ — > | and |^ — > | contribute almost equally into total. If one 
considers the full contribution of the interference terms, namely from the |^ — > | and 
1^ — > I transitions, one get just about 30% of the total squared terms. Note, that in 
total interference terms contribute destructively. To within error bar accuracy, one may 
conclude that the dB{El)/dE distribution for ^^C is composed of an incoherent sum of 
several squared terms that peak at about the same energy, corresponding, to the 2si/2 <S> 2"*" 
component of the g.s. configuration. This attests to the validity of the simple cluster model 
even in this subtle quantal system. 

In Ref. Q, the final channel was clearly identified as ^^C(2''")+n. Therefore, to be 
correct, one has to compare the data with the exclusive Coulomb dissociation cross-section, 
Eq. (fTK|) . However, as we have seen, the core 0+ state contribution is at most 10% of the 
peak value of the 2"*" contribution. Further, it is mostly concentrated in the relative energy 
range < E < 2 MeV. Thus one may safely calculate the inclusive cross-section and ignore 
the low energy part. The comparison with the data is shown in Fig|3 It can be seen, that 
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FIG. 3: Contributions to the dipole strength function for C. (a)All squared contributions to 
the dipole strength for |^ — *■ | ; (b)Dominant interference contributions for |^ — > | ; (c)All 



contributions to the dipole strength for 4 



3 + 



^ ; (d) Comparison of overall contribution of 



squared terms (contribution of O''' is not shown here) and interference terms. 



the present calculation reproduces the data. 
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FIG. 4: C+ Pb at 495A MeV Coulomb breakup cross sections calculated in different models. 



Data from 



E* is an excitation energy. (Here the theoretical cross sections were convoluted 



with the experimental response functions.) 
C. ^^C 

The last example is the calculation of the Coulomb breakup of ^^C. Strictly speaking, we 
cannot use the rotor model for the core, as it was experimentally shown, that the first excited 



22] • But, as it was shown in Ref. 



state of the core is not 2+ but 1~ at around 7 MeV 
particle-rotor model reproduces the spectroscopic number rather well. This calculation gives 
87% of the core inert state. Because of this the contribution of core excited components is 
very small and we do not include it here, that is we only include contributions 2si/2 — > P3/2 
and — > Pi/2 to the dipole strength function. The resulting cross section is shown in 
FiglHl along with the result of the cluster model Eq. lfTSj) . Again the simple cluster picture 
works nicely for this system. 



V. CONCLUSIONS 



We have used the particle-deformed core model wave functions to calculate the Coulomb 
dissociation of ^^C, ^''C and ^^C in the semiclassical model. We have shown that whereas 
interference effects are important in the non-halo isotope ^''C, they are practically insignifi- 
cant in ^^C and ^^C. This fact points to the validity of the cluster picture in the latter cases 
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FIG. 5: C+ Pb at 605A MeV Coulomb breakup cross sections calculated in different models. 



Data from 



E* is an excitation energy. (Here the theoretical cross sections were convoluted 



with the experimental response functions.) 

where it emphasizes the subtle quantal nature of the former, mostly core-excited projectile. 
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